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Although ion channels are attractive targets for drug discovery, the systematic screening of ion 
channel-targeted drugs remains challenging. To facilitate automated single ion-channel recordings for the 
analysis of drug interactions with the intra- and extracellular domain, we have developed a parallel recording 
methodology using artificial cell membranes. The use of stable lipid bilayer formation in droplet chamber 
arrays facilitated automated, parallel, single-channel recording from reconstituted native and mutated ion 
channels. Using this system, several types of ion channels, including mutated forms, were characterised by 
determining the protein orientation. In addition, we provide evidence that both intra- and extracellular 
amyloid-beta fragments directly inhibit the channel open probability of the hBK channel. This automated 
methodology provides a high-throughput drug screening system for the targeting of ion channels and a 
data-intensive analysis technique for studying ion channel gating mechanisms. 

Membrane proteins, such as transmembrane receptors and ion channels, are major molecular targets for 
commercial drug therapies Compared to the successes achieved in the efficacies of drugs targeting 
transmembrane receptors, such as G protein-coupled receptors, only 30% of the drug-targeted ion 
channels predicted by in silico gene expression analyses are being utilised in current therapies^. When searching 
for drugs that target ion channels, it is important to identify the site of drug action and to determine whether the 
site is intra- or extracellular. This type of topological information facilitates the determination of the duration of 
drug administration, thereby reducing side effects^'^. Moreover, single- channel recording of point-mutated ion 
channels is essential to determine how each drug affects channel conformation. However, topological (intra- and 
extracellular) single ion-channel recording is labour intensive and requires highly skilled researchers; even then, it 
is difficult to acquire an extensive data set that can be used to reveal biologically relevant molecular mechanisms. 
Therefore, an automated parallel measurement technique for single-channel recording is required for efficient 
drug discovery or ion-channel analysis. 

In general, two electrophysiological methodologies can be used to examine the effect of a ligand on the gating 
mechanism of an ion channel: first, the conventional inside-out and outside-out patch-clamp method using living 
cells^'^, and second, the artificial bilayer lipid membrane (BLM) method^"^^. Single -channel recordings have not 
been achieved for either method using an automated system. For the patch-clamp method, an auto-patch system 
has been commercialised but does not provide inside-out/outside-out recordings to identify the intra- and 
extracellular effects of drugs at the single-molecule leveP^. In contrast, the BLM method of ion channel analysis 
provides a platform for achieving the equilibrium recognition of the intra- or extracellular topology using purified 
materials at the single-channel leveP^. Although several attempts to create the parallel formation of BLMs have 
been reported^^"^\ single- channel recording for identifying membrane topology has not been achieved in the 
parallel system due to the low stability and reproducibility of BLMs. Therefore, these conventional methods 
present challenges in the preparation of a reproducible bilayer and its stability for use in automated, parallel 
recording^^. 

Here, we report an automated, parallel methodology for single ion-channel recording based on Droplet 
Contact Method (DCM)^^. In this method, a BLM is formed at the interface between contacting droplets sub- 
merged in an oil/lipid mixture (Fig. la, b), which is recently referred to as Droplet Interface Bilayer (DIB)^^. To 
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Figure 1 | Automated, parallel ion channel-recording system using the droplet contact method, (a) Formation of a bilayer lipid membrane by the 
droplet- contact method. Two lipid monolayers are brought into contact with one another and spontaneously form a bilayer. (b) The PDW chip for the 
simple preparation of lipid bilayers. (c) Schematic illustration of the PDW chip, (d) Schematic diagram of the PDW chip connected to a patch-clamp 
amplifier, (e) An arrayed 16-ch PDW chip connected to a multi-patch- clamp amplifier, (f) Automated, parallel bilayer single- channel recording system. 
The injection robot automatically prepares the bilayer lipid membranes, reconstitutes the ion channels and controls the drug application. 



enhance the stability of the BLM, we reduced the area of the BLM by 
placing a polymer sheet with microsized pores in between the drop- 
lets (Fig. Ic). This improvement allows us to form highly stable and 
reproducible BLMs. Our system has advantages over the previous 
system in simultaneously satisfying the abilities as following: i) 
intensive data acquisition, ii) single- channel analysis, and iii) analysis 
of topologically identified targets. Using this method, we performed 
single- channel recordings on several types of ion channels with auto- 
mation for intensive data acquisition, and we identified whether 
these channels had intra- or extracellular sites of drug action. We 
also demonstrated that amyloid-beta fragments (AP) directly inhi- 
bits the opening of a Ca^^ -dependent (hBK) channel at both the 
intra- and extracellular domains. 

Results 

Automated lipid bilayer preparation using an injection robot. A 

double-well chip (Fig. Ib-d) was used for the ion channel recordings 
in conjunction with the droplet contact method as previously 
reported^^ In this approach, lipid/oil and buffer solutions are 
sequentially deposited onto a double-well chip to prepare the BLM; 
this method has been broadly used as a simple means of BLM 
formation^^'^^. Here, the stability of the lipid membrane on the 
chip was improved for use in the automated system (as discussed 
below). To decrease the membrane capacitance and to increase its 
mechanical stability, the contact area of the lipid monolayer was 
reduced using a poly(chloro-p-xylylene) (parylene) separator that 
was patterned with micropores^^'^^. The parylene film was 
integrated into the double-well chip to produce a parylene double- 
well (PDW) chip (Fig. Ic). Therefore, the BLMs that formed in the 
PDW chip were mechanically stable as described in the following 
section. Due to the high stability of the BLMs, we were able to 
construct parallel BLMs using an injection robot; the resulting 
arrayed BLM device, which consists of a 16-channel (16-ch) PDW 
chip, is shown in Figure le and f and Supplementary Movie L 

Stability and noise examination of BLMs prepared by the droplet 
contact method. We first estimated the mechanical stability of the 



BLMs. The mechanical stability of a BLM prepared by the traditional 
method cannot be tested because it is too fragile^^'^^. However, due to 
the robustness of the BLMs formed with the DCM, mechanical 
stability was evaluated using a mechanical oscillator. The survival 
rate of the BLMs in the double-well (DW) chip lacking the parylene 
separator and the PDW chip as a function of the vibration frequency 
is depicted in Figure 2a. The survival rate of the DW chip was less 
than 40% at a vibration of 30 Hz. In contrast, the BLMs in the PDW 
chip maintained their form within this frequency range. The 
resonant frequency of the BLMs in the PDW chip was calculated 
to be 1.3 kHz (Supplementary Text 6). In the range of ambient 
vibration frequencies encountered in the daily environment 
(<100 Hz)^°, the BLMs in the PDW chip were mechanically stable. 

Electrical current noise is an important aspect of artificial cell 
membrane experiments^\ Figure 2b shows the power spectra of 
the mean current noise of the ot-haemolysin (otHL) open channel 
signal in the DW chip and the PDW chip. Due to the small mem- 
brane capacitance of the PDW chip, the current noise at lower fre- 
quencies (<100 Hz) was two orders of magnitude lower for the 
PDW chip than for the DW chip. The overall membrane capacitance 
was reduced with the reduction of the bilayer area, which resulted in 
the reduction in the current noise. 

The long-term stability of the BLMs on the PDW chip was also 
measured by obtaining alamethicin channel signals (Fig. 2c). Since 
this stability is required for long-term experiments such as cell-free 
expression experiment, the evaluation of the long-term stability is 
important. The signal event frequency was approximately 10 Hz in 
this experiment (Supplementary Text 6). Although the event fre- 
quency was slightly reduced compared with the initial state, the 
channel signals were observed for 14 days; compared to previous 
reports, this period is the longest reported for BLM stability^^'^l 

Introducing ion channels into the planar lipid bilayer using a vesi- 
cle fusion technique. Vesicle fusion is an efficient methodology for 
channel reconstitution in planar lipid bilayers^^. The single- channel 
current of the purified and reconstituted ion channel on the surface 
of a proteoliposome can be observed using this method. The process 
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Figure 2 | (a) Mechanical stability in the DW and PDW chips as 
determined with a mechanical vibrator. The survival rate was defined 
as the number of stable BLMs/16-ch for 30 s = 4, mean ± S.D.). 
(b) Comparison of the power density spectra of aHL open-channel 
currents between the DW and PDW chips. The PDW chip showed 
significantly reduced noise in the lipid bilayer. (c) The long-term stability 
of BLMs in the PDW chip. Channel opening events are presented as a 
function of days. The channel signals of alamethicin were observed for 14 
days. 

of vesicle fusion is described in Supplementary Text 5. The fusion 
rate of the vesicle containing the hBK channel was defined as So (So = 
signal-obtaining chamber/ 16-ch), and it was used to optimise vesicle 
fusion in the PDW array (n = 16,N = 3). We found that So increased 
with increasing hydrostatic pressure (Supplementary Text 7) at the 
interface between the droplets (Fig. 3a). This result can be explained 
because the area of the bilayer formed in the droplet interface is 
increased by the high pressure, thereby increasing the probability 
of vesicle fusion. Meanwhile, in the higher pressure region, the 
BLMs were ruptured within 120 min. The optimal pressure was 
determined to be 14 N m~\ So was stable at 0.6 under these 
conditions after 120 min. 
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Figure 3 | The optimisation of vesicle fusion, (a) Comparison of the hBK 
channel signal rate (signals/ 16-ch) at the first signal observation and 
after 120 min with various hydrostatic pressures at the droplet interface 
(AT = 3, n = 16, mean ± S.D.). (b) The signal rate was increased with 
increasing concentrations of proteoliposomes. The single-channel signals 
overlapped and could not be discriminated at levels higher than the 
optimal concentration (AT = 3, /7 = 16, mean ± S.D.). 

The other important factor for vesicle fusion is the concentration 
of the vesicle. Figure 3b shows the concentration dependence of S^. 
The curve relative to the concentration was convex upward. At levels 
under 10"^ mg/mL lipid, the probability of fusion between the vesicle 
and the planar bilayer should be enhanced by increasing the vesicle 
concentration. In contrast. So was decreased at levels higher than the 
apical concentration. The current signals overlapped due to an excess 
of vesicles that fused after the initial vesicle fusion. Therefore, So was 
decreased in the high-concentration range. The highest So value 
(more than 0.9) was obtained at a lipid concentration of 10"^ mg/ 
mL. In addition, the addition of cholesterol to the BLM decreased So 
(Supplementary Text 7). 

Previously, few studies of vesicle fusion probability have been 
made on the basis of quantitative data due to the requirement of 
labour-intensive experiments. Here, the vesicle fusion conditions 
were optimised based on extensive data that were obtained using 
the automated BLM method. The optimised condition of the channel 
reconstitution was utilised in the following experiments. 

Verification of the automated system for measuring several ion 
channels and a mutant channel. To test the device, toxin channels 
(aHL and alamethicin), channels (hBK and Wild-Type 
(WT)KcsA) and a mutant channel (KcsA mutant (E71A)) were 
reconstituted in the 16-ch PDW chip using the automated system 
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Figure 4 | Examples of single- channel current recordings using the automated system, (a-e) Current traces (left) and corresponding point 
histograms (right) from (a) a-hemolysin (holding potential (Vhoid) = +50 mV), (b) alamethicin (Vhoid = +150 mV), (c) hBK (Vhoid = +60 mV), 
(d) WT KcsA and (e) E71A ( V^oid = + 100 mV). The arrows indicate the channel closed state, (f, g) /-^curves of the hBK channel (f) and KcsA channels 
(WT and E71A, g). The data are shown as the mean ± S.D. (n > 3 each). 



(Supplementary Text 5). E71 is an amino acid residue of the pore 
helix region of KcsA, and its corresponding mutant (denoted E71A 
here) has been used in studies on channel-gating mechanisms^^. The 
currents that passed through each ion channel were observed in 
parallel at the single -channel level under various recording 
conditions. Figures 4a-e show the current signals and histograms 
of five different proteins at the single- channel level. The l-Y curves 
from the hBK, WT-KcsA and E71A channels were measured from a 
single run of the automated system (Figs. 4f and g). Compared to the 
WT-KcsA channel, E71A displayed an increased channel open 
probability (Po) (Figs. 4d and e) and a decreased outward 
rectification due to a reduced energy barrier for ion entry into the 
selective filter at negative holding potentials'^ These differences 
between the WT and the mutant are consistent with previous 
reports^^. The application of this methodology to gene-targeted ion 
channels will provide essential information regarding drug-targeted 
amino acid residues that contribute to the molecular basis of drug 
action at specific sites on ion channels. 

Massive data acquisition for the hBK channel using automated 
and parallel single-channel recordings. Parallel recordings of the 
hBK channel were performed using the 16-ch PDW chip. The hBK 
channel controls firing patterns in neurons and modulates the tone 
of blood vessels^^. The current recordings from the PDW chip are 
shown in Figure 5a. Rapid open-and-close transitions associated 
with channel gating were observed at the single-channel level from 
all of the channels in the PDW chip. Because this parallel 
measurement easily facilitated a high degree of reproducibility, we 
obtained numerous single- channel recordings {n = 102) and used 
the large amount of data to analyse the channel conductance and 



dwell time (Figs. 5b and c). The large amount of data obtained 
ensured reliability in the measurement of these biophysical pheno- 
mena. Furthermore, a particular channel orientation could be obser- 
ved by controlling the solution contents. For the hBK channels, the 
current signals were observed from oriented channels, in which the 
Ca^^ -binding site was directed toward a Ca^^ solution that was added 
to one of the compartments of the PDW chip. The characteristics of 
the oriented hBK channels were investigated, and the results are 
shown in Supplementary Figure S4. The fundamental properties of 
the hBK channels, which were measured with the automated parallel 
system, were similar to previously reported results obtained with live 
cells and BLMs'^'^^ 

Topologically identified hBK channel inhibition using the auto- 
mated system. The greatest advantage of the automated system is 
that it can be used to examine the direct inhibition of purified 
channels from both membrane topological directions (intracellular 
and extracellular) while excluding interactions with cytosolic com- 
ponents. Here, we demonstrate the ability to analyse direct inter- 
actions between the ion channel and ligands that are applied from 
either the intra- or extracellular side. Iberiotoxin (IbTX) is a toxin for 
K^ channels and inhibits the hBK channel from only the extracellular 
side by decreasing Po^°- This extracellular inhibition was observed at 
the single- channel level in our system (Supplementary Fig. S6c, d. 
Supplementary Text 9). The dose-dependent inhibition of IbTX on 
Po was evaluated to determine the IC50 and the Hill coefficient, which 
were 5.6 nM and 2.1, respectively (Fig. 6a). In addition, the inhi- 
bition of hBK by the L-type calcium channel blocker verapamil 
was also examined. The dissociation constant (k^) for verapamil 
was determined to be 6.3 jiM (Fig. 6b, Supplementary Text 9). 
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Figure 5 | The performance of the 16-ch recording system, (a) Parallel recordings of the hBK channel obtained using the automated system. 
Single- channel currents were observed from all 16 channels of the arrayed PDW chip. The arrows indicate the channel closed state, (b, c) Histograms 
of the conductance and dwell time of hBK channels obtained using numerous single- channel recordings (n = 102, each recording period > 20 s). 
The holding potential was +60 mV. 



These measurements were consistent with previous reports on the 
inhibitory effects of IbTX and verapamil on reconstituted hBK 
channels^°'^\ This type of pharmacological analysis using our auto- 
mated system has great potential for aiding in the determination of 
drug activity. 

Amyloid-beta fragments directly inhibited a single hBK channel. 

We investigated the topological effects of an amyloid-beta fragment 
(Api.42) on hBK at the single-channel level. Figures 7a and c show the 
current recordings obtained from the hBK channel before and after 
the injection of APi.42 on the intracellular side. The value 
decreased immediately after the injection. A similar reduction was 
also observed following an extracellular application (Fig. 7d). The 
effect of A p 1.42 on the lipid bilayer itself was also examined because 
amyloid proteins can form a channel structure in negatively charged 
membranes (DOPE)^^. However, APi.42 did not affect the neutral 
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Figure 6 | Automated single hBK channel inhibition from both the intra- 
and extracellular sides, (a) Dose-dependent inhibition of the hBK channel 
by iberiotoxin (IbTX) and (b) verapamil (n = 4, mean ± S.D.). 



membrane used in this study (EggPC) as shown in Figure 7b. This 
experiment demonstrates another advantage of this artificial system: 
the lipid content can be controlled to directly observe channel 
activities without disturbing the membrane. Figures 7e and f show 
the current histograms that were obtained before and after the 
application of A Pi. 42 to a single hBK channel from the intrace- 
llular (Fig. 7e) and extracellular (Fig. 7f) sides. In both cases, the 
open channel states were inhibited by APi.42. The effects of APi.42 
on channel conductance are shown in Figure 7g; the same buffer 
solution without A p 1.42 was used as a control (Ap (-)). The con- 
ductance of the hBK channel did not change upon A Pi. 42 appli- 
cation. In contrast, Pq decreased upon A Pi. 42 application from 
both directions (Fig. 7h). 

Discussion 

Our automated methodology represents a powerful means to acquire 
a massive amount of data regarding topologically identified single 
ion channels. Using this method, we investigated hBK channel 
inhibition with respect to membrane topology (intracellular or extra- 
cellular) by using the amyloid-beta fragment (Api.42) . Abnormal BK 
channel activity has been recently reported in the pyramidal cells of 
Alzheimer's disease model mice^^. However, the effect of AP on the 
hBK channel remains unclear for two reasons. First, it is unknown 
whether AP affects the BK channel directly or via another in vivo 
inhibitory cascade. Second, it is unknown whether Ap affects the 
channel from the intra- or extracellular direction. Although extra- 
cellular deposits of A P in the cortical neuropil and blood vessel walls 
in the brain are classical histopathological features of Alzheimer's 
disease, emerging evidence from transgenic mice and human 
patients has revealed that Ap also accumulates intracellularly in 
neurons^^. Therefore, we used our system to measure the inhibitory 
response of Api.42 applied to topologically identified hBK channels. 
It was revealed that A Pi. 42 directly inhibits the hBK channel from 
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Figure 7 | Automated single hBK channel inhibition by APi.42 from both 
the intra- and extracellular sides, (a) Effect of Api.42 on single hBK 
channels before and after the addition of A Pi. 42 to the extracellular side, 
(b) Effect of Api.42 application on the lipid bilayer (EggPC) itself, (c, d) The 
intra- (c) and extracellular side (d) effects of the appUcation of APi.42 to the 
hBK channel, (e, f) Point current histograms for the hBK channel obtained 
before (top) and after (bottom) the application of APi.42 from the 
intracellular (e) or extracellular side (f) direction, (g, h) Effects of Api.42 on 
the channel conductance (g) and open probability (h) (n = 12, mean ± 
S.D.). The significance was assessed by t-tests: *P < 0.05, **P < 0.01. 

both the intra- and extracellular sides. These results suggest that in 
neurons, A Pi. 42 directly inhibits the hBK channel at both the intra- 
and extracellular domains. A similar inactivation by A Pi. 42 may also 
occur with other ion channels, such as extracellular AP-sensitive ion 
channels^^'^^; this process may represent one cause of neuronal dys- 
function in the brain of Alzheimer's disease patients. 

Taken together, our automation system readily revealed the 
topological and direct effects of various compounds, including 
pharmacological agents and APi.42, on ion channels at the single- 
molecule level. In the DCM system, the inner solution of the droplet 
can also be exchanged to perform drug application assays^^. 



Moreover, automated ion-channel screening using the droplet con- 
tact method may also be used to analyse mutated ion channels, in 
which the determination of the site of drug action is necessary for 
personal gene therapy. Further comprehensive studies with this 
automated ion channel system will contribute to the systematic dis- 
covery of drugs that target ion channels. 

Methods 

Automated lipid bilayer preparation by the droplet contact method and ion 
channel introduction into BLM. The following method was used to form the bilayer 
lipid membrane (BLM) by the droplet contact method: 1) An EggPC/ ;i-decane 
(20 mg/mL) solution (6 to 9 uL) was initially deposited into each well. 2) The aqueous 
buffer and the same buffer containing the water-soluble membrane proteins or 
proteoliposomes (15 to 20 uL) were then deposited into each well (Supplementary 
Fig. S2a). The lipid monolayers contact each other and form a BLM within 10 minutes 
of the addition of the buffer solution. The BLMs can be formed in the array device 
using an automated injection robot (epMotion 5070, Eppendorf). A movie of the 
automated preparation of BLMs is provided (Movie SI). 

We defined the two droplets that sat on the recording and the ground electrodes as 
droplets A and B, respectively. A symmetrical buffer solution (in mM: 1000 KCl, 2 
KH2PO4, 8 K2HPO4, 1 EDTA, pH 7.4) was used in both droplets A and B for the 
water-soluble toxin channels (ocHL and alamethicin). ocHL and alamethicin were 
dissolved in droplet A at 600 and 10 nM, respectively. In contrast, the reconstitution 
of the hBK channels and KcsAs was more difficult because the channels are water- 
insoluble and retain their conformation only within a lipid bilayer membrane. Hence, 
the reconstitution of these channels was conducted using vesicle fusion. The process 
of vesicle fusion involves the following steps: i) A vesicle (the ion channel- expressing 
proteoliposome) approaches a BLM (Supplementary Fig. S4a). ii) The vesicle and the 
lipid bilayer fuse with each other when the vesicle attaches to the bilayer surface. To 
enhance the fusion event, the osmotic pressure between the inside and outside of the 
vesicle was controlled in this experiment (Supplementary Fig. S4c). iii) The channel 
signals are measured through Ag/AgCl electrodes after the reconstitution. The 
experimental conditions used for the simultaneous and parallel recordings are 
described in Supplementary Text 5. In the binding experiments, a small volume of 
solution containing the inhibitors can be injected into the droplet after determining 
the single-channel currents. This type of kinetic experiment was impossible using 
previous microdrop methods. 

Ion channel current recordings and data analysis. The channel current signals were 
detected using a multi-patch-clamp amplifier (JET-Bilayer, Tecella) with a 1 kHz 
low-pass filter at a sampling frequency of 5 kHz (unless otherwise noted) at 23 ± 1°C. 
Most of the experiments on the ion channels, including the pharmacological tests, 
were performed using a K-recording buffer (in mM: 140 KCl, 0.1 CaCl2, 10 HEPES, 
pH 7.0) or an Na-recording buffer (in mM: 120 NaCl, 0.1 CaCl2, 10 HEPES, pH 7.0). 
The current analysis was performed using the pCLAMP ver. 10.6 (Molecular Devices, 
Sunnyvale) and Igor Pro 6.2 (Wavemetrics, Oregon) software programs. The single- 
channel analysis of the ion channels is described below. The channel recordings were 
conducted for 2 h after injecting the proteoliposome. If the BLM was ruptured, it 
could be reformed by repairing the BLM using a hydrophobic needle. The channel 
signals typically began approximately 1 to 10 min after injecting the droplets. At least 
20 s of continuous recording was used to estimate the channel open probability (Pq). 
The durations of open and closed times were measured with a half- amplitude 
threshold analysis using the pCLAMP software program. Occasional long closed 
periods (shut intervals > 5 s) were excluded because of extremely low activity. To 
promote a decreased signal/noise ratio, an off-line digitised low-pass filter (<500 Hz) 
and the baseline offset were used. The results are given as the mean ± S.D. 

In case of the topology identification of the channels, the channel which is oriented 
undesired direction can be discriminated with changing the voltage polarization 
because of the rectifying. For instance, using voltage- dependent KcsA channel (as 
shown in Fig. 4), the oppositely oriented channel shows the inappropriate conduc- 
tance under applying the opposite voltage. As the result, we can identif)^ the channel 
orientation by the channel conductance. On the other hand, if the channel does not 
have the rectifying, the topology can be recognized by observing the direction of the 
blocker or the opener application. In addition, the specific orientation can be main- 
tained owing to the conical shape of the transmembrane domain, which matches the 
high curvature of small vesicle (<500 nm dia.)"**^. The orientation is discriminated 
from the direction of the vesicle fusion. 
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